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In this article we present first an algorithm for calculating the determining equations associ- 
ated with so-called "nonclassical method " of symmetry reductions (a la Bluman and Cole) for 
systems of partial differential equations. This algorithm requires significantly less computation 
time than that standardly used, and avoids many of the difficulties commonly encountered. The 
Q I proof of correctness of the algorithm is a simple application of the theory of Grobner bases. 

C/2 ■ In the second part we demonstrate some algorithms which may be used to analyse, and 

often to solve, the resulting systems of overdetermined nonlinear pdes. We take as our principal 
, . example a generalised Boussinesq equation, which arises in shallow water theory. Although 

rS I the equation appears to be non-integrable, we obtain an exact "two-soliton" solution from a 

j^ ■ nonclassical reduction. 
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1 Introduction 

Nonlinear phenomena have many important apphcations in several aspects of physics as well as 
other natural and applied sciences. Essentially all the fundamental equations of physics are nonlin- 
ear and, in general, such nonlinear equations are often very difficult to solve explicitly. Consequently 
perturbation, asymptotic and numerical methods are often used, with much success, to obtain ap- 
proximate solutions of these equations; however, there is also much current interest in obtaining 
exact analytical solutions of nonlinear equations. Symmetry group techniques provide one method 
for obtaining such solutions of partial differential equations (pdes). These have many mathemat- 
ical and physical applications, and usually are obtained either by seeking a solution in a special 
form or, more generally, by exploiting symmetries of the equation. This provides a method for 
obtaining exact and special solutions of a given equation in terms of solutions of lower dimensional 
equations, in particular, ordinary differential equations (odes). Furthermore they do not depend 
upon whether or not the equation is "integrable" (in any sense of the word). 

Symmetry groups have several different applications in the context of nonlinear differential 
equations (for further details see, for example, [S, p^] and the references therein): 

• Derive new solutions from old solutions. Applying the symmetry group of a differential 
equation to a known solution yields a family of new solutions (quite often interesting solutions 
can be obtained from trivial ones). 

• Integration of ODES. Symmetry groups of ODES can be used to reduce the order of the 
equation (such as to reduce a second order equation to first order). 

• Reductions of pdes. Symmetry groups of pdes are used to reduce the total number of 
dependent and independent variables (for example, reduce a pde with two independent and 
one dependent variable to an ode). 

• Classification of equations. Symmetry groups can be used to classify differential equations 
into equivalence classes. 

• Asymptotics of solutions of pdes. Since solutions of pdes asymptotically tend to solutions 
of lower-dimensional equations obtained by symmetry reduction, some of these special solu- 
tions will illustrate important physical phenomena. Furthermore exact solutions arising from 
symmetry methods can often be effectively used to study properties such as asymptotics and 
"blow-up" . 

• Numerical methods and testing computer coding. Symmetry groups and exact solutions of 
physically relevant PDES are used in the design, testing and evaluation of numerical algorithms; 
these solutions provide an important practical check on the accuracy and reliability of such 
integrators. 

The classical method for finding symmetry reductions of pdes is the Lie group method of in- 
finitesimal transformations and the associated determining equations are an overdetermined, linear 
system (cf., [^, |3l|]). Though this method is entirely algorithmic, it often involves a large amount of 
tedious algebra and auxiliary calculations which can become virtually unmanageable if attempted 
manually, and so symbolic manipulation programs have been developed to facilitate the calcula- 
tions. An excellent survey of the different packages presently available and a discussion of their 
strengths and applications is given by Hereman [p^. 
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There have been several generalizations of the classical Lie group method for symmetry reduc- 
tions. Ovsiannikov p^] developed the method of partially invariant solutions. Bluman and Cole |M], 
in their study of symmetry reductions of the linear heat equation, proposed the so-called "nonclas- 
sical method of group- invariant solutions" (in the sequel referred to as the nonclassical method), 



which is also known as the "method of conditional symmetries" |25| and the "method of partial 



symmetries of the first type" [44|. This method involves considerably more algebra and associated 



calculations than the classical Lie method. In fact, it has been suggested that for some pdes, the 
calculation of these nonclassical reductions might be too difficult to do explicitly, especially if at- 
tempted manually since the associated determining equations are now an overdetermined, nonlinear 
system. Furthermore, it is known that for some equations such as the Korteweg-de Vries equation, 
the nonclassical method does not yield any additional symmetry reductions to those obtained using 
the classical Lie method, though there are pdes which possess symmetry reductions that are not 



obtained using the classical Lie group method. Olver and Rosenau [B3, 34] proposed an extension 
of the nonclassical method and concluded that "the unifying theme behind finding special solutions 
of PDES is not, as is commonly supposed, group theory, but rather the more analytic subject of 
overdetermined systems of pdes" . 

Clarkson and Kruskal [^] developed an algorithmic and direct method for finding symmetry 
reductions (hereafter referred to as the direct method) and using it obtained previously unknown 
symmetry reductions of the Boussinesq equation. The novel characteristic about the direct method 
in comparison to those mentioned above, is that it involves no use of group theory; additionally, 
for many equations the method appears to be simpler to implement than either the classical or 
nonclassical methods. 



Olver ||3j (see also [@, |36|]) has recently shown that the direct method is equivalent to the 
nonclassical method when the infinitesimals for the independent variables are autonomous with 
respect to the dependent variables, generating a group of ^''fibre-preserving transformations" . 

It has been known for several years that there do exist pdes which possess symmetry reductions 



that are not obtained using the classical Lie group method (cf., |3^, |3^). Recently the direct 



and nonclassical methods have been used to generate many new symmetry reductions and exact 



solutions for several physically significant pdes, which represents important progress (cf., |lj, 21] 
and the references therein). 

In §2 we present an algorithm for calculating the determining equations for the nonclassical 
method for a system of pde. These are usually calculated by reducing the so-called infinitesimal 
equations, obtained from the group prolongation of E, with respect to both S and the invariant 
surface conditions ^. However, in practise, reducing an equation with respect to a system is not 
well-defined. Indeed, unless one is careful about the choice of term in each equation from which 
to back-substitute, infinite loops can occur in the reduction process. The theory used to overcome 
these difficulties is that of Grobner bases, a powerful computational tool in algebra, geometry and 
logic. Our essential idea is to first reduce the given system T, using the invariant surface conditions 
^ to generate a simplified system S^. Then we apply the classical Lie method to S^. We use the 
theory of Grobner bases, as they apply to algebraic systems, to provide a proof of correctness of an 
algorithm for finding the determining equations for the nonclassical method which eliminates the 
problems and moreover proceeds efficiently. 

In §2.1 we describe the process of applying the nonclassical method, and the difficulties en- 
countered, in more detail. Then we give an elementary introduction to Grobner bases adapted 
to our purposes, and show why they solve the problems. Next we give the algorithm for finding 
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the determining equations and prove it is correct, followed by some worked examples that are 
prototypical. 

Since "nonclassical symmetries" of S are actually classical symmetries of the system consisting 
of T, augmented by the invariant surface conditions, the basic idea of using Grobner bases to make 
the reduction of the infinitesimal equations well-defined applies to finding classical symmetries of 
any system. We stress that the "nonclassical symmetries" obtained by the nonclassical method are 
not symmetries of S itself, since they do not necessarily transform all solutions of the system to 
other solutions. Rather, they are symmetries of the augmented system given by S together with 
specified auxiliary conditions. 

In §3 we discuss some algorithms and strategies that have proved useful in making solving 
the determining equations for the nonclassical method, which are overdetermined and nonlinear, 
tractable. We illustrate this with some examples. 

The difficulties of the nonclassical method due to the determining equations being an overde- 
termined nonlinear system makes the use of symbolic manipulation programs more important. 
Levi and Winternitz |25| and Clarkson and Winternitz [18| in their applications of the nonclassical 



method to the Boussinesq and Kadomtsev-Petviashvili equations respectively, interactively used 
the MACSYMA program SYMMGRP.MAX [^]. Nucci |3C] has also developed an interactive program 



NUSY in REDUCE for the nonclassical method. Here, we use the maple package diffgrob2 |26(| , 
which appears to be the only differential algebra package available that can handle equations not 
solvable for their leading derivative term. In the appendix we give details of how our algorithm to 
obtain the determining equations may be implemented using symmgrp.max. The interesting thing 
here is that we are using the SYMMGRP.MAX program for a purpose for which is was not originally 
designed, however, it can be adapted to generate the determining equations for the nonclassical 
method since the latter can be interpreted as the determining equations for the classical method 
applied to an appropriate system of equations. 



2 Determining Equations for Symmetries 

2.1 The Classical and Nonclassical Methods 

Suppose one is given a system of partial differential equations 

S = {/i = 0,...,/, = 0}, (2.1) 

where each /j is some polynomial expression in the independent variables {xi, . . . ,Xn}, the depen- 
dent variables {ui, . . . , Um} and the derivative terms {uk^a I ol G N"} where 



dx"^' . . . dx^- ' 



Uk,a = TTal ET-^- (2-2) 



One also can have transcendental and arbitrary functions of the dependent variables in what follows 
without affecting the theory. 

We recall briefly the method of finding the determining equations for classical Lie point symme- 
tries thereby fixing our notation. Let u^ ' denote the list Ufc^o-, where k = 1,. . . ,m and \a\ = N. 
The index a + i is given by (ai, ...,«» + 1, ... , «„), while ct + 7 = (ai + 71, ...,«„ + 7n)- 
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To find the classical Lie point symmetries of the system S, one takes a group action defined 
infinitesimally by 

X* =Xi + eCi{x,u) +0{e'^), i = 1, 2, . . . ,n, 
u* =Uj + e4)j{x,u) + 0{e'^), j = 1,2, . . . ,m, 

where x = xi, . . . , x„ and u = ui, . . . , u^- Then one requires that this transformation leaves the 
set of solutions 

5s = Maj)|/i = 0,/2 = 0,...,/, = 0}, 

invariant. 

The N^^ order partial derivatives transform according to (with |q;[ = A^) 

f)N * 
a *,ai *,a2^ STJTT =Uj^a+e(t)^f^ (x,u,u^^\ . . . ,u^^A +0{e^), 

OX{ X2 . . .Xn -^ ^ ^ 

where the A^*'^ extension, denoted by (/>'• [x,u, u^^', . . . , u^^' j , is given by the recursive formula 



,[a] 



=1 



■^ \ ) \jxi r^ Dx,- 



and 






is the total derivative operator [31|. We make the obvious definition 






^Q: — j^ ai y.^a2 ■ ■ ■ j^^a„ ■ (2-5) 



Now consider the system 

{f,[x*,u*,u*^'\x*),...,u*^''\x*))=0\i = l,2,...,r], 
which is ( p.ip with u replaced by u* and a; by a;*. It is easily seen that 

fi(^x*,u*,u*^'\x*),...,u*^''\x*)) 

= f,(x,u,uW{x),...,u^''\x))+epr^^^^^{f,) + 0{e'), 

where 

j=l J fc=l '^ fc=l |q;|>1 ' 

is known as the N^^ prolongation (or A^*'^ extension) of the infinitesimal operator 



^ = YCjix,u)^— + YMx,u) — 

j = l ^^0 fc=l '^"'^ 
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Let A*" be the order of the system S. Requiring that (^]^) is invariant under the transformation, 
i.e. 

pr(N)v(/,)^_ =0, i = l,...,r (2.7) 

yields an overdeter mined, hnear system of equations for the infinitesimals ^{x,u) and 4>{x,u), 
obtained by setting the coefficients of the different monomials in the derivative terms {uk,a \ 1 ^ 
k < m, \a\ 7^ 0} in the pr^^-* v(/j) = to zero. This means that the N^^ prolongation of /j, 

i = l,2,...,r, is zero whenever u is a solution of the original system (^.1[). 



The important point to note is that we are considering the invariance under the group action 
of the system viewed as an algebraic system in the indeterminants {xi,Uj,Uj^a}, in the relevant jet 
bundle (|31|, Ch. 2). 

The nonclassical method of Bluman and Cole Q for finding symmetry reductions of a system 
of PDES involves appending the invariant surface conditions to the system and finding the classical 
symmetries of the appended system. The invariant surface conditions are given by 

-fpi = ClUi^xi + C2Ui,X2 + ■ ■ ■ + (.nUi,x„ - 4>i = 0, i = l,2, ...,m. (^) 

In this method one requires only the subset of Sj] given by 

S^^ = {u{x, t)|/i = 0, /2 = 0, . . . , /, = 0, Vi = 0, 7/^2 = 0, . . . , ^^ = 0} , 

is invariant under the infinitesimal transformation. The idea is that since the invariant surface 
conditions map to themselves under the prolonged group action, they are not a restriction on the 
infinitesimal equations of the system S, but rather since ui^xi and wi,x2) for example, are no longer 
independent, the determining equations will be more general than those for the classical method. 
Hence imposing these conditions leads to the possibility that there are more solutions, not less. 

The usual approach to finding the determining equations for the nonclassical method of the 
system S consists of calculating the infinitesimal equations pr'^^ v(/) for / G S, where A^ is the 
order of the system S, and reducing them with respect to the augmented system G = S U ^. 
By reduction with respect to G is meant elimination (or back-substitution) from pr^'^) v(/) of 
derivatives of one pre-determined term from each of the equations in G. One then reads off the 
coefficients of the different monomials in the derivatives of the Uj; setting these to zero are the 
determining equations. 

In practise there are several difficulties. 

Example 2.1.1. Consider the equation 

Ai = uxt - f{u) = 0. (2.8) 

The infinitesimal equation pr'^' v (Ai) (with xi = x, 3:2 = t, ^ = ^i, ^2 = 1 and <j) = 4>i), is 

(jjM _ f'(u)(t) = (pxt + 4>xuUt + {4>tu - ixt)Ux + {(t>uu - Cxu)UxUt 

-ituui - S.uuUtul + {(pu - ix)Uxt - itUxx - '^iuUxUxt - CuUtUxx " f'{u)(j). 

When reducing this equation before reading off the determining equations, does one reduce the 
derivative term Uxt using the original equation, or using the x derivative of the invariant surface 
condition 

ilj = i{x,t,u)ux + ut-(l) = Ql (2.9) 
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The difference of the two reductions is proportional to a differential consequence of the system 
{Q, Q}, namely, 

—-^p - Ai = f{u) + {^x + CuUx) Ux + ^Uxx - (Px- 4>uUx = 0. 

Dx 

Using this equation, one can eliminate all Uxx terms, given ^ ^ 0. But should one? This leads to the 
next question: "By which differential consequences, if any, do we need to reduce the infinitesimal 
equations in order to obtain the determining equations for the nonclassical method?" 

Example 2.1.2. Second, consider the equation. 

As = utt - Uxx = 0. (2.10) 

With the same notation as the previous example, pr^^' v (A2) is: 

(Pxx ~r y^Yxu t^xxj'U'x ~r \(puu ^<,xu)Ur[ s,uuU^ + \(pu ^'ixj'Uxx "Jsu^x^xx 

-{4>U + 20t„nt - ittUx + 4>uuUt - '^ituUtUx - iuuU^Ux + (l)uUu - "^itUtx - iuUttUx - 2^uUtxUt}- 

Instinctively one would eliminate the t derivatives using the invariant surface condition ut = (f)—^Ux, 
and then eliminate the Uxx terms using Uxx = utt- But eliminating Uxx introduces a utt term, and 
eliminating this leads to a Uxt term, and eliminating that leads to a Uxx term again. Clearly 
care must be taken in the reduction procedure (by which is meant the successive eliminations or 
back-substitutions) to prevent infinite loops occurring. 

These difficulties are all part of the problem of finding the "normal form" of an equation / 
with respect to some given system of equations G. By "normal form" is meant some well-defined 
reduction of / such that no further eliminations from G are possible. 

For questions concerning "normal forms" to be answered in a well-defined way, the concept of a 
Grobner basis is required. The algorithm to calculate a Grobner basis for a system of polynomials 
over a field was developed by Buchberger p, ^J and since then has been extended to a wide variety 
of algebraic scenarios; the concept has a large number of applications (see [IC, 19] and references 



therein). In the following section, we discuss Grobner bases as they are used in our particular 
application. 

2.2 Grobner Bases for Differential Polynomials 

Consider the set 

3 = {xi,ei,^,0j-^,Uj-a I 1 < i < n,l < j < m,5 G N"+™,a G N", !«!, \5\ < N}, 

where N is some finite number. The system of equations S, and the infinitesimal equations ob- 
tained by the prolongation of the group action on S, can be considered to be polynomials in the 
elements of 3 with complex coefficients. In this section, we discuss Grobner bases for systems of 
polynomials. We give those definitions required for the sections that follow, and examples, relevant 
to our application. An excellent introduction to Grobner bases can be found in ||n 



We denote the set of polynomials in the indeterminates {Ci, . . . Xs} with coefficients in C by 
TZ = C[Ci, ■ ■ ■ ,Cs]i and use the multi- index notation for multiplication defined for (3 = (/3i, . . . , /5s) G 
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N'^ by C = Ci^ • • • Cs" ■ Take an ordering O on the indeterminates, say Ci < C2 < • • • < Cs, and 
define the lexicographic ordering lex(C?) on the set of monomials {C^ \ (3 G N*} to be 

C^ >lex(0) C' (2.11) 

if for some < j < s — 1, 

13s = 5s, ■■■, (3s-j+i = 5s-j+i, 13s- j > 5s- j- (2-12) 

Following Bayer and Stillman [H], another monomial ordering denoted here by BS(r) is given by 

^ >BS(r) '' ^^ ^^ + ^''+1 "^ V I3s> 5r + (5r+i H h 5^, else for some j, [3s = Ss, • • • , /3s-j+i = 

(5s -j 4-1, /3s -j < <5s-j- There is a wide variety of monomial orderings available |^]. They are 
required to have the so-called compatibility, or multiplicative, property, that is, 

C/^>C^^>CC'^>CC^ and C'^C^>C^. (2.13) 

Given an ordering on the indeterminates {Ci}' ^^ have defined above several orderings on the 
monomials {C }• We now discuss orderings on the indeterminates 3- 

We assume that Uj^a > (Pk 6 ^ ^i6 ^ ^i ^^^ ^^^ •?' ^' ^' *' ^ ^^^ ^^ ^^^ require that the ordering 
chosen on our indeterminates is compatible with the differential structure, that is, we take the 
ordering on the indeterminates to be such that 

Uj,a > Uk,f3 =^ Uj,a+j > %,/3+7 and Uj,a.+'y > Uj^a, h\ / (2.14) 

and similarly for the ^ and the 4>. 

In Example 2.1.2 where the reduction process was infinite, an incompatible ordering had been 
chosen. In that example we were eliminating t-derivatives using S^Ux + ut — 4> = 0, which implies 



ut > Ux- Then by compatibility (2.14), we must have that uu > Uxt and u^t > Uxx, so that 



utt > Uxx- Thus using Uxx — utt = to eliminate occurrences of Uxx is an incompatible choice. 

Finally, we require that derivatives of the u with respect to one particular pre-chosen variable, 
Xk (say), are greater in the order than derivatives with respect to other independent variables. This 
is needed to apply the elimination ideals property of Grobner bases |^^ in the proof of correctness 
of our algorithm. Thus we want a compatible ordering which is decided first with respect to 
some ordering on the dependent variables, say Um > Um-i > • • • > ui, and then the number of 
derivatives with respect to x^, and then any compatible choice thereafter. Such an ordering on the 
indeterminates we will denote by the term /c-order, or 0{k). For a k-order 0{k) choose r = r{k) 
such that for I > r, Q represents a derivative term Uj^a with ak 7^ 0, and for I < r, Q represents 
either a derivative term Uj^cx with a^ = or one of the remaining indeterminates. We then denote 
the monomial ordering BS(r) by BS(C(/c)). 

Example 2.2.1. For second order systems of PDE in two independent variables, suppose we 
want a compatible ordering on the derivative terms of u such that ut > Ux- Two possibilities are 
utt > Uxt > Uxx > Ut > Ux > u and utt > Uxt > ut > Uxx > Ux > u. The first is an ordering where 
total degree of differentiation is one of the deciding factors in the ordering, the second ordering is 
a t-order. 

To find the coefficient of a monomial M in p, denoted coef(p, M), one looks for all summands 
in p of the form rjM where tj E C. Then one defines coef(p, M) = J2 n- H coef(p, M) 7^ 0, we say 
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that M occurs in p. The highest monomial term occurring in a polynomial p is denoted hmt(p) 
and its coefficient is denoted hcoef(p). 

Example 2.2.2. Consider the polynomial 

P = iu4>xUul^UxxUt + {ixx - iu)UttUxxUx + x'^{<f>u - S,x)UxtUtt + U^xtUxxUttUt ■ 

The independent variables are x, t, the dependent variable is u, the level of prolongation is 2, so that 
0{x) is given initially by Uxx > Uxt > Ux > utt > ut, and 0{t) is given initially by uu > Uxt > ut > 
Uxx > Ux- The orders 0{k) are not unique. Then in the ordering lex(C(x)), iiMT{p)=uttUxxUx^xx, 
in lex(C'(t)) hmt {p)=u%UxxUt^u(t>xU while in BS(0(t)), iiMT{p)=UxxUuu^uCxt- 

Definition 2.2.3. Suppose for some polynomial q that the HMT(q) divides some monomial M that 
occurs in the polynomial p, so that ^ HMT(g) = M, and that hcoei{q) = a £ C. Then the reduction 
of p at M with respect to q is given by 

p^,p-coe{if,M)C^q/a. (2.15) 



Thus reduction depends on the ordering used. The use of a compatible ordering ( p. 13 , p. 14 ) 



eliminates the infinite loops observed possible in the Introduction; see P, 27] where it is proved that 



with respect to a compatible ordering, reduction is a noetherian relation, that is, it must terminate 
after a finite number of steps. 

The definition of reduction uses no differentiation, since we wish to remain within the algebraic 
domain. Suppose A^ is the highest degree of differentiation occurring in the given system S. When 
considering the system S U ^, the ^ equations need to be prolonged to order N, and the system 
to which our theory will apply will be 

G = S U {Bai; = | V G 'I', « G N", \a\ < N - 1} (2.16) 



where Dq. is defined in (2.5). We assume that all equations in G are of the same order of differen- 



tiation; if not, those of lesser order need to be prolonged. 

Definition 2.2.4. A normal form of a polynomial p with respect to a system of polynomials G, 
denoted normal(p, G), is achieved when no further reduction of p with respect to any member of G 
is possible. 

Definition 2.2.5. The ideal generated by a finite system of polynomials G C C[Ci, . . . , Cs], is the 
set I{G) = {EgeGfg9\ fa ^C[Ci,...,Cs]}. 

Definition 2.2.6. A Grohner basis of an ideal I{G) is a finite set H C C[Ci,...,Cs] such that 
I{H) = I{G) and where for all p G I{G), one has normal(p, H) = 0. Thus a Grobner basis depends 
upon the ordering used. We denote a Grobner basis for the ideal generated by G in the monomial 
ordering order to be Gb(G, order). 

Sufficient conditions to characterise a GB and an algorithm to calculate the GB for an ideal over 
a field in any compatible ordering were given by Buchberger p], and it and its generalisations can 
be found now in textbooks (see for example |]l^]). This algorithm has since been implemented in 



various symbolic algebra programs, for example in mathematiga [^, maple ||13|], reduge ||2 
or the specialist package magaulay ^. 

Given a Grobner basis G of a polynomial ideal I{G), the normal form of any polynomial with 
respect to G is well-defined. Let / be a polynomial, and let hi and /i2 be two normal forms of / 
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with respect to G. Then /ii — /12 € I{G) by the reduction formula. Since G is a Grobner basis for 
I{G), hi — /i2 reduces to zero with respect to G. Since neither hi nor h2 reduces with respect to 
G, neither can their difference (reduction must take place at some monomial which must occur in 
at least one of hi or in /i2), so that difference must already be zero. Thus different reductions of / 
with respect to a Grobner basis are equal. 

Example 2.2.7. For the standard example 

Uxx = A{x,t,u,Ux,ut,Uxt,utt), (2.17) 

with invariant surface condition 

Cux + ut = (^, (2.18) 

one needs to prolong the invariant surface condition to be of the same order as the given equation, 
namely 2, so one has the system G: 

Uxx- Mx,t,U,Ux,Ut,Uxt,Uu) =0, 

(Cx + iuUx)Ux + iUxx + Uxt - (t>x - (f>uUx = 0, 

(6 + iuUt)Ux + iUxt + Utt - 4>t- (t>uUt = 0, 

S,Ux+Ut- (j) = 0. 

Take the ordering given initially by uu > Uxt > Uxx > ut > Ux- Finding the Grobner basis of these 
equations is equivalent to considering the first three of them as equations in Uxx, Uxt and uu and 
calculating the "echelon" form of the system. In the ordering given initially by uu > Uxt > ut > 
Uxx > Ux, we eliminate all the ut terms using ^Ux + ut — t] = from the other conditions, and then 
calculate the echelon form of the reduced system to obtain conditions for uu , Uxt and Uxx ■ 

To prove the correctness of our algorithm to generate the determining equations for the non- 
classical method, we need the elimination ideals property of Grobner bases ||l^, ^ : 

Theorem 2.1 Let C[Ci, . . . , Cr] be the set of all polynomials over C in the first r indeterminates. 
Suppose G is a Grobner basis of the ideal I{G) in the ordering lex(O) where O is Ci < C2 < 
. . . < Cs \2.11 , 2.1^ . Then for for all 1 < r < s, G D C[Ci, . . . ,Cr] generates the elimination ideal 



/(G)nc[Ci,...,Cr]. 

In words, this property means that for every r, any polynomial in the first r indeterminates that 
can be found from the generators by addition and by multiplication by elements of C[Ci, . . . , Cs]) 
can be "read off' from the Grobner basis. 

For the Bayer and Stillman ordering BS(r), with r pre-deter mined, we have that GnC[Ci, . . . , Cr] 
generates /(G) n C[Ci, . . . Xr], that is, we obtain only the one elimination ideal |5|, |l^. For our 
application, we only need one particular elimination ideal, that given by using BS(0(A;)), while the 
Bayer and Stillman orderings are more efficient than the lexicographic orderings j^. 

2.3 The Algorithm to calculate the determining equations for the nonclassical 
method. 

Let S be a system of pde of "polynomial type", with n independent variables {xi, . . . , Xn}, and m 
dependent variables {ui, . . . , Um}- Terms like sm{uj) or h{uj) where h is undetermined or arbitrary 
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are simply considered to be additional indeterminates in the algebraic computations and do not 
affect what follows. 

There are n cases to consider, namely where the invariant surface conditions {^) have the forms, 
for 1 < j < m, 

^iUj,xi+C2Uj,x2 + ■ ■■ + Cn~iUj,x„.i+Uj,xn =4>j-, Case n 

6'"i,xi + 6%,X2 + ■ • • + %,a:„_i =(t>j-, Case n - 1 



'^j x-[ — T*? ; v^ase 1. 

That is, we have successively for 1 < /c < n, that ^^ has been set equal to 1 and ^k+i = • • • = ^n = 0. 
We consider each case separately. In the following algorithm, in the A:*^ case, we require 0{k) to 
be a fc-order, and order to be one of lex(C'(A;)) or BS(0(A;)). 

Algorithm: Determining equations for the Nonclassical Method 

INPUT: S, a system of pde with highest order of differentiation A^; 

A; G {1, . . . , n}; 0{k), a /c-order, orderg {lex(0(fe)), BS(0(A;))}. 
OUTPUT: DetSqns, the determining equations for the nonclassical method of the system S 

in the k^^ case 

for j from 1 to m let 

ipj := ^lUj^xi + • • • + 6-i^j,Xfe_i + Uj^xk - 4>j 
^* ■= {Da^|Jj I 1 < j < m, Q € N", |q;| < A^ - 1} 
/C:={normal(/,^*)|/ES} 
Tn/:={pr(N)v(/)I/e/C} 
QB :=gb(/C, order) 
ninf := {normal(/, gi3)\f e Inf} 

VetSqns := {coef(/, ui,a, . . . Um,(xJ = | / G TZInf, otj G N" \ {0}} 
end 

In words, we reduce the given system S with respect to the invariant surface conditions, form the 
infinitesimal equations of the result, reduce the infinitesimal equations with respect to an algebraic 
Grobner basis of the reduced system, and then read off the coefficients of the result in the usual 
manner. 

By reducing the equations in E with respect to \&*, that is, eliminating all derivatives of the 
Uj with respect to xt-, before finding the infinitesimal equations, the Grobner basis calculation is 
greatly diminished. Indeed, for systems consisting of a single equation, the calculation is elimi- 
nated altogether, and in that case the algorithm becomes the classical method, but on the reduced 



equation. Since the algorithm to calculate a Grobner basis has poor complexity |25], our method 
is more efficient. 

Before discussing the correctness of the algorithm, let us demonstate it on Example 2.2.8, 



Example 2.3.1. Given the general second order equation ( 2.17] ), with invariant surface condition 
(|2.18 ), the algorithm splits into two cases, r ^ and r = 0. 



Case I. T ^ 0. In this case we set r = 1 (without loss of generality), and eliminate ut^ Uxt and uu 
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in (|2T^ ) using (l2l8|) , i.e. 



ut 

Uxt 

utt 



h - iu^ (2.19) 

^x + ^uUx - {ixUx + ^u'"^ + iu,j,,j) (2.20) 

^t + 0«(</' - C'Wx) - 6'"x - CuUxifj) - S.Ux) - ^[(Px + <t>uUx - {ixUx + iuU^. + ^Ma;a;)]-(2.21) 



Substituting these into ( 2.17|) yields the ode (with t a parameter) 

A(x, t, u, Ux,Uxx;^, ix-,it-, £,u, <P, <t>x, 4>t, 4>u) = 0. 



(2.22) 



Now apply the classical Lie algorithm to this equation. Thus we apply the second prolongation 
pr(2)v to ( 2.22| ) and require that the resulting expressions vanish foru€5 = ■^u:A = 0> =5^, 



I.e., 



pr 



(2) 



v(A) 



A=0 



0, 



where ^ and (f) appear both in ( ^.22 ) and pr'^V. Equating coefficients of powers of Ux to zero then 
generates the determining equations. 

Case II. r = In this case we set ^ = 1 and so the invariant surface condition reduces to 
^x = (pixjtju). Hence we obtain the differential consequences 

Uxt = (t>t + 4>uUt, Uxx = 4>x + (t>uUx = (t)x + 4>4>x- 

Substituting these into (|2.171 ) yields the ode (with x a parameter) 

A(x, t, u, Ux,ut; 4>, (px, (pt, (pu) = 0. (2.23) 

Now apply the classical Lie algorithm to this equation. Thus we apply the second prolongation 
pr^^V to ( 2.23|) and require that the resulting expressions vanish for u G 5 = <u{x,t) : A = Ok 
i.e., 



pr 



('K(a 



A=0 



0. 



Equating coefficients of powers of ut to zero then generates the determining equations. 

Proof of Correctness: We apply the theory described in §2.1. We have that S is a system of 
PDE in the dependent variables {ui, . . . , Um} and the independent variables {xi, . . . , Xn} of order 
A^. In the /c*^ case, we have 



n 

and G 



{6%,xi + i2Uj,x2 + ■■■ + Uj,xk - '/'j I 1 < J < m}, 

{DaV-lV e ^,a eN", \a\ < iV - 1}, 

C[xk,Uj,(t).^f^,^^^f^Uj,a\ l<k<n,l<j<m,l3e N"+'", a G N", |q|, |/3| < N], 

s u ^* c 7^. 



By definition, in the notation used in this article, the determining equations for the nonclassical 
method of the system S are 

jcoef (normal(prW v(/), Gb(S U ^*)),ui,ai ■ ■ ■ Urn,a^) = | / G S, «, G N" \ {0}| (2.24) 



The following Lemma is proved in [|31| (Proposition 3.33 and Theorem 3.38). 
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Lemma 2.3.2 The group prolongations of the invariant surface conditions satisfy pr v(^) E /(^*) 
for ip G I(^*). 

Lemma 2.3.3 In the A;*^ case of the algorithm, where ^* = {Dcxip \^p £ "^,0^ e N", |a:| < iV — 1} 
and /C = {normal{f, ^*) [ / € S}, then pr*-^-* ^rg, for g £ IC have all derivative terms of the form 
Uj,a satisfying Ofc = 0. 

la] 
Proof: Consider the term (f)- . By examining the formulae for the prolongation of the group 



action, (2[^,2.6), there are two ways that a derivative of Uj can occur in the pr^^ ^((7): if either a 
derivative of ^^ is non-zero, or if g has a term of the form Uj^a. with a^ 7^ 0. Since ^^ is constant, 
and since no equation in S after reduction by ^* contains a derivative of Uj with respect to Xk (for 
any j), the lemma is proved, n 

Let / G S and let /' denote normal(/, ^*). By Lemma 2.3.2 and the formula for reduction, 
noting that we have hcoef(-(/') = 1 for all V' G 'I'*, /' = / + Z^^'e** ^4''^ implying 

pr(N) v(/') = prW v(/) + K (2.25) 

where r^ & TZ and k G /(^*). Hence reducing both sides of ( |2.25| ) with respect to Gb(I; U ^*) 
yields the same result. However, in reducing / by ^*, the choices of the r^ have removed the need 
to reduce the left hand side of ( |2.25| ) by ^* at all. This is the content of Lemma 2.3.3. Hence there 
is no need to reduce the elements of Inf by the invariant surface conditions. 

Li *&* we have a set of equations for the Uj^a^k occurring linearly. Now a GB obtained from a 
reduced set of generators is still a GB so that gb(/C, order) is equal to the elimination ideal 

GB(G, order) n C[uj^a I |a| < ^, a/c = 0, 1 < i < m]. 

Hence it is sufficient to reduce the elements of In/ with respect to Gb(/C, order). □ 

Since the system is regarded as an algebraic system on the relevant jet bundle, we see why 
differential consequences obtainable only by further differentiation, the so-called integrability con- 
ditions, are not relevant. Of course, one can always investigate the result of reducing with respect to 
additional integrability conditions [37|. The following remark by Olver and Rosenau |34] is relevant 



here: ". . . the reason why Bluman and Cole find nontrivial conditions on their groups in order to 
apply their nonclassical method is that they fail to take into account the additional restrictions on 
the derivatives of u coming from . . . integrability conditions." According to ^^, every group is a 
"weak symmetry group" if all the integrability conditions are taken into account. 

Another idea is find the integrability conditions first and then calculate the determining equa- 
tions for the enlarged system, discussed by Schwarz |42| for classical symmetries. 



2.4 Examples 

In this section we give some prototypical examples of what our method looks like in action. The 
two examples cover those cases mentioned as being problematic in the Introduction. In the Ap- 
pendix, we give the input files with which the macsyma package symmgrp.max [|12| calculates 
the determining equations for the second example. Provided at some point a program written 
to calculate classical symmetries recognises that the terms representing infinitesimals in the input 
equations are the same functions used by it to represent the infinitesimals, the correct equations will 
be obtained. For this it may be necessary to input the infinitesimals in the internal representation 
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used for them by the package. Since we are using the package for a purpose for which is was not 
originally designed, it is important to know precisely how it can be adapted. 

In this section we set xi = x, X2 = t, ^i = ^ and ^2 = t- 

Example 2.4.1. The Nonlinear Wave Equation Uxt — f{u) = 0. 

For the nonlinear wave equation in characteristic coordinates 

Uxt = fiu), (2.26) 

we use the differential consequence of the invariant surface condition with r = 1 to eliminate Uxt 
in ( 2.26| ) and so obtain 

0x + (puUx - {ixUx + iuul + iuxx) = f{u). (2.27) 

Now apply the second prolongation pr'^^V to this equation and eliminate Uxx using (2.27). Actually, 



one can only eliminate ^Uxx] if one's program eliminates Uxx so that the equation being reduced 
is multiplied by ^, that is equivalent to putting ^ in the coefficient ring and to assuming that ^ is 
non-zero. In some cases it may be necessary to consider the case ^ = separately. In this case, 
allowing ^ = is equivalent to the next case, r = 0. Finally, equating coefficients of powers of Ux 
to zero yields the following four determining equations: 

iiuix - i'^ixu - itiu - 4>uiiu + ^itu + 4>uui^ = 0, 

f{u)iix - 4'4'xCu - f{u)(t)^u - 4>xit - f{u)it + 4>4>xuC - f{u)(puC + <Pxt^ + f'{u)(l)i = 0, 

<t^iuix + itix - (t^iixu + 4>xiiu + 3/(n)C^„ - 4>(t)uiu - iixt - (puS.t - (pxuC"^ + 4>4>uuC + 4>tu( 



0. 



In the case r = we set ^ = 1 and so the invariant surface condition reduces to Ux = (j){x, t, u). 
Then we use the differential consequence of this to eliminate Uxt in ( 2.26| ) and so obtain 



f{u) = (f)t + (t>uUt- 



(2.28) 



Now apply the first prolongation pr^^V to this equation and then eliminate ut using ( p. 28 ). This 
yields one determining equation 



/) + <t>XU {(t>t - f) - (t>U<t>: 



•xt 



iJuVtu 



+ 



5n3- = 0. 

au 



Example 2.4.2. A generalised Boussinesq equation. 

Here we derive the determining equations for a "generalised" Boussinesq equation 

Utt + Uxx + aUxUxt + PutUxx + Uxxxx = 0, (2.29) 

where a and /3 are arbitrary nonzero constants. This equation, together with several variants, can 



be derived from the classical shallow water theory in the so-called Boussinesq approximation |46]. 



Furthermore the Painleve PDE test due to Weiss, Tabor and Carnevale |45], suggests that the 
equation is non-integrable for any non-zero choice of a and /?. 



In the case when r ^ 0, we set r = 1. Using the invariant surface condition ^Ux + ut 
eliminate ut, Uxt and uu in ( 2.29| ) to yield 



we 



h + 4>u{4>- S,Ux 

+ Uxx + aux 



■ itux - iuUx{(t> - iux) -i (i)x + 4>uUx - {ixUx + iuui + iUx 

+ (t>uUx - {ixUx + ^uul + iUxx) + PUxx (</» - (Ux) + Uxxxx 



0, (2.30) 
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Apply the classical Lie algorithm to this equation using the fourth prolongation pr^^W and eliminate 
Uxxxx using (2.30), to yield the following determining equations, 



(puu = 0, 
^(pxu ^Kxx ^ U, 

-2l3(,(j)xu - ai(t>xu + cifku + aix4>u + a4>tu + I^Hxx - 2a^^ - a^xt = 0, 
-a^cAx + Q<i>xxu + P(p(pu + /?</>* + 2p^x^ - A^xxx + 4.fCx + 2ix + 2iit = 0, 

(t>xxxx + P(p4>xx + (t>xx + acpl - 4:^Cx(t>x " 2^4 (/'x + ^ix4>(t>u + 4>U + '^(fxptu + 4^x0* = 0, 
PC(l>xx - '2^0i(l)u4>x - Oi^x(t)x - ^xxxu - '^.P^Xpxu - a(l)(j)xu " 2</>a;w + ^iix(t>u 

+2^t4'u - a(t>xt + 2i4>tu + Kxx(t> + ixxxx + ixx - ml + 2itix + itt = 0. 

In the Appendix we give the macsyma input files used to calculate these determining equations 
using the package symmgrp.max. We discuss the solution of these determining equations and the 
associated symmetry reductions in Example 3.2.2 below. 

In the case when r = 0, we set ^ = 1. Using the invariant surface condition Ux = <j){x,t,u) we 
eliminate Ux, Uxx, Uxt and Uxxxx in ( ^.291) to yield 

utt + a(p{(j)t + 0«nt) + (1 + (3ut){(j)x + (PM 

+ (f'xxx + 4'u4>xx + 3(t)(t)x(t>uu + ^(t>x(t>xu + 4>x(t>l + 4>^(t>uuu 

+ '^(f^'Pxuu + 4^(t>'^(t>u(f>uu + 3(f>4>xxu + 5(t>4'u(t>xu + (fxt't = 0. (2-31) 

Applying the second prolongation pr^^' v to this equation yields the following determining equations; 

^uu = 0, 

P(l)xx + (« + /?) {(px^u + Hi) +{a + 2/3)00^„ + 204„ = 0, 

(pxxxx + 4>xx + (t>U + 2(j)ct)xu + ^xu(t>xx + Hx(t>xxu + 4:(j)4>uxxx + ^x4'u4>xu + Q<t>4>u4>xxu 

+A(pfyxu + 8(j)(pl^ + (q + /3) {4>x(pt + (I^Mu) + acp {(pxt + (p(ptu) = 0. 



3 Algorithms for solving the systems of determining equations 

3.1 Introduction 

Finding the determining equations for the nonclassical method is only half the story. Since the sys- 
tems are overdetermined and nonlinear, it is necessary to have additional algorithms and strategies 
to aid in their analysis and solution. In this section we demonstrate the use of the Kolchin-Ritt 
algorithm |11, ^ |2^, ^, ^] in conjunction with the DirectSearch [|^] and Reid |4C] strategies to 



solve the systems of determining equations. 

In §2 of this article, algebraic procedures were used. The process of calculating algebraic 
Grobner bases of successive prolongations of a system (see for example [|ll[) for the purpose of 
finding integrability or compatibility conditions is impractical, even if one has an a priori bound on 
the level of prolongation necessary to find all such conditions. Instead, we employ a true differential 
analogue of Buchberger's algorithm, with cross-multiplication replaced by cross-differentiation, and 
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algebraic reduction replaced by differential reduction, and so on. For fully nonlinear systems, that 
is, containing equations that are not linear in their highest derivative terms, it is necessary to 
use pseudo-reduction instead of reduction else the differential algorithm will not terminate. The 
resulting algorithm is called the Kolchin-Ritt algorithm and has been implemented in packages in 
MAPLE (H, 0). 

In the Kolchin-Ritt algorithm, each pair of equations in the system given or obtained en route is 
cross-differentiated so that their highest derivative terms become equal. One then cross-multiplies 
by the relevant coefficients so that these terms cancel. The result is then pseudo-reduced by all 
known conditions, and if non-zero, is called an integrability or compatibility condition of the system. 
This process continues until no new conditions are found. In pseudo-reduction one may multiply 
the equation being reduced by nontrivial terms in order to effect the reduction. A simple example 
will illustrate our meaning. Suppose one wishes to reduce the equation / : Uxyy — UxxUy = by the 
equation g : UyUxy — u^ = 0. In lex(C'(x)), we have that the highest derivative term in g is Uxy 
Then a (one-step) pseudo-reduction of / with respect to g is given by 

In strict reduction, one is allowed to multiply the equation being reduced by at most expressions 
in the independent variables. The use of pseudo-reduction means that the build up of differential 
coefficients needs to be taken into account when interpreting the results of the calculation. Precise 
statements describing the output and the limitations of the Kolchin-Ritt algorithm can be found 
in 1^^, with additional algorithms that eliminate some (but not all) of the problems. If the output 
equations of Kolchin-Ritt are linear in their highest derivative terms, then one will have obtained 
either the integrability conditions or the differential analogue of the elimination ideals of the system, 
depending on whether the ordering used is total degree or lexicographic. In other cases, more work 
may be required to obtain the maximum amount of information possible. 

The Reid strategy is designed to overcome as far as possible the problem of the build up of 
differential coefficients by adjusting the order in which pairs are chosen to be cross-differentiated. In 
this strategy, the system is divided into four groups of increasing complexity; single term equations, 
linear equations, equations linear in their highest derivative terms, and fully nonlinear equations. 
Into the fourth class is also put those equations with excessively many terms or high degree of 
differentiation. Each class is reduced with respect to those in the lower classes. Beginning with 
the linear class, the compatibility or integrability conditions are calculated and the result used to 
(pseudo)-reduce the equations in the higher classes. This can lead to nonlinear equations becoming 
linear for example. If a condition is found that belongs to a lower class, the strategy is to recalculate 
the algorithm on the lower class with the new condition added. If a condition is found that belongs to 
a higher class, the strategy is to put it there and "keep it for later". In this way the need to multiply 
by differential coefficients in the reduction process is delayed for as long as possible, and kept to an 
absolute minimum. This strategy has been implemented in maple in the StandardForm procedure 
of the Reid-Wittkopf Differential Algebra Package ||4^. We give an example to demonstrate this 
strategy. Example 3.3.2. During the course of the calculations, if any equation is of the form of an 
expression raised to a power, we throw away the power to keep the complexity of the calculation 
down. In addition, we throw away factors known to be non-zero, but we do not throw away factors 
involving arbitrary parameters. 

The DirectSearch algorithm is a strategy that searches specifically for conditions less than 



those in a given set in the ordering in use. This strategy was used effectively in |16] to solve a 
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classification problem. The first example demonstrates the use of the DirectSearch strategy to 
make the calculation of the Kolchin-Ritt algorithm tractable. 

3.2 Examples 

We start with a simple example, namely the determining equations for the nonclassical method 
for Burgers' Equation. We then study a more complicated example, the generalised Boussinesq 
equation of Example 2.4.2. All the calculations in this section were performed using the maple 
package diffgrob2 ||2^. By a lexicographic ordering based on ui < U2 . . . < Um and xi < X2 < 
. . . < Xn is meant an ordering on the derivative terms such that Uj^ct < '^k 3 '^^ ^ ^ ^ ^^^^ ^^ ^ (^^ 
else an-i < Pn-i and so on. It is assumed that Uj^a > x^ (of course). Arbitrary constants of 
integration in this section are denoted by Kj. 

Example 3.2.1. The determining equations for the nonclassical method in the r = 1 case for 
Burgers' equation 

(3.1) 

are 



20^« 



Calculating the Kolchin-Ritt algorithm on this system leads to expression swell, so we use a little 
finesse. Take a lexicographic ordering based on </> > ^ and x > t > u. In this ordering, the 
determining equations are listed in ascending order. The DirectSearch on the first three equations 
consists of the following procedure. We cross-differentiate the second and third equations so that 
the leading terms cancel, and then reduce with respect to ^uu = 0. We iterate this procedure 
recursively using the result of the previous calculation and the second determining equation until 
we have eliminated <J3. We then continue the process using the first determining equation until we 
get zero. This procedure terminates yielding 

Cn(Cn + l)(2Cn-l)=0. 

Thus not only is S,u a constant, we know that it can take at most three values. Taking each value 
of ^u in turn yields three tractable calculations. 

Case 1. Reducing the determining equations with respect to ^^ = and performing the Kolchin- 
Ritt algorithm yields, after simplification, 

^u = 0, 

(/) + <„ + 2^ex + 6 = 0. 

The equations have been listed in ascending order. The aim of using a lexicographic ordering is 
to try to reduce the integration of the system to that of a series of ode by starting with the least 



Ut = Uxx + '^UUx 




SM« ~ 


= 0, 


Ku - 2Cux + 2eC« + 2uCu -- 


= 0, 


'2(f>ux + ixx - "^ixi - ixu - it = 


= 0, 


(t^xx - 4>t- '^.(pix + U(t>x = 


= 0. 
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equation and successively integrating up. Setting ^ = G{t)x + H(t) into the third equation and 
equating coefficients of powers of x to zero yields 

G"{t) + 6G{t)G'{t) + AG^t) = 0, 
H"{t) + 2G'{t)H{t) + 4G^{t)H{t) + AG{t)H'{t) = 0. 

The first can be linearized to z"'{t) = by setting G{t) = z'{t)/2z{t) while the second is linear in 
H given G so that 



G{t) = ^ /^ + ""^ ^ H{t) = I K3 

Thus we have i^, and given S, the fourth equation gives 



8(2t + Ki) 4 



+ KAG{t). 



Case 2. Reducing the original determining equations with respect to ^„ + 1 = and performing 
the Kolchin-Ritt algorithm yields = and ^ = —u, in other words the solution is output! 

Case 3. Reducing the second of the determining equations by 2^„ — 1 = we obtain (puu+C+u = 0) 
implying 2(/)ut(u +3 = 0. Thus (/> is cubic in m. Substituting = —^u^ + G{x,t)u^+H{x,t)u+K{x,t) 
into the remaining equations, eliminating S, using (/>„„ + ^ + u = 0, and setting the coefficients of 
different powers of u to zero yields the system S for G, H and K, 

—Gt + Gxx + '^GxG + Hx = 0, 

Kxx-Kt + AKGx = 0, (3.2) 

Hxx + Kx-Ht + AHGx = 0. 

In a total degree ordering, the Kolchin-Ritt algorithm on this system yields no new integrability 
conditions, and the formal general solution space appears to depend on six arbitrary functions. Since 
the system is highly symmetric in its derivative terms, is nonlinear and is not overdetermined, we 
need to proceed with caution. Reverting to a lexicographic ordering, with K > H > G and t > x, 
cross-differentiating the second and third equations and reducing yields 

AKGxx + '^Hxxt — Htt + SHtGx + AHGxt — Hxxxx — 8HxxGx — SHxGxx — ^HGxxx — IQHG^ = 0. 
Thus, unless Gxx = 0, we have solved for K. We first solve this singular case by calculating the 



Kolchin-Ritt algorithm on (3^) together with Gxx = to yield. 



'^xx — U, 

-Gut + AGxuG + 12GxtGt + l2GxGu-^8GxGxtG-48GlGt + MGlG = 0, 

Hx-Gt + 4GxG = 0, 

-Htt + SHtGx + 4HGxt - WHGl + 2Gxtt - 24G,G,t + 32Gl = 0, 

Kx-Ht + AHGx + Gxt-4Gl = 0, 

-Kt + AKGx + Gtt-4GxtG-8GxGt + l6GlG = 0. 



(3.3) 



Setting G{x, t) = xa(t) + b{t) in the second equation of ( |3.3| ) and reading off coefficients of x yields, 

120^2 - attt - 96ata'^ + lea^a + 64a^ = 0, 
4aff 6 - bttt + I2atbt + 126^0 - 486*0^ - 48atab + 640^6 = 0. 

The first of these can be linearised to C""{'t) = by setting a{t) = —('(t)/{4:({t)), and substituting 
this into the second equation yields (6(t)C(i))'" = 0. Using this and the final four equations in the 
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output (^), which are hnear in H and K, yields the singular solution. This solution generalises 
that in Q]. 

We return to the general case of (|3.2|), and use the condition obtained for K above to eliminate 
K. Using a Direct Search method we can solve for H in terms of G, and then obtain two very 
long expressions for G whose cross-derivative is zero. None of the singular cases that arise lead to 
solutions other than that obtained in the Gxx = case. This is as far as elementary differential 
algebra can take one with (|3.2| ), and further analysis will depend on additional input, such as any 
geometrical information that comes from knowledge of what the G, H and K actually mean. This 
could lead to a transformation in which the analysis is simpler and the geometric content clearer. 
One possibility is the connection, for Burgers' equation, between the nonclassical method and the 
so-called Singular Manifold method of obtaining solutions using a truncated Painleve expansion 



discussed in [20|. This uses an ansatz in which G satisfies a Burgers' type equation, so that the 
process of finding the nonclassical reductions becomes a Backlund transformation for the Burgers' 
equation. Another ansatz, used in j^, is Gj = 0. However, the full solution appears to be unknown 
at present. Further, the meaning of such a large solution space of nonclassical reductions appears 
to be unknown. 

Example 3.2.2. In this example we apply the classical and nonclassical methods for symme- 
try reductions of a "generalised Boussinesq equation" ( 2.29| ). The classical symmetry reductions 



(3.4) 



are calculated not only for comparison purposes, but to show just how effective the Kolchin-Ritt 
algorithm can be for linear systems. We do not consider the case where either a and P is zero. 

3. 2. 2. A Classical Method. Applying the classical Lie point method to this equation yields the 
system of linear overdetermined equations given below (with ^i = 6. and ^2 = t)- 

Tu = 0, (ptu - ixt = 0, 

Tx = 0, Tt- 2ix = 0, 

Cu = 0, 'Pu-Tt + 2Cx = 0, 

4>uu = 0, (2/? + a)(j)xu - Pixx = 0, 

2</>a;« - "iixx = 0, fi(t)xx + '2(t>tu - Tu = 0, 

(/? + a)it = 0, (l>xxxx + (Pxx + <t>u = 0, 

a(l)x - 2it = 0, Q<i)xxu + P^t - ^ixxx + 2ix = 0, 

^(Pxxxu ~r ^(pxu ~r OL(pxt C,xxxx <,xx stt — U. 

Before applying the Kolchin-Ritt algorithm to a linear system, the set is simplified to be an "auto- 
reduced" set. In the diffgrob2 package the procedure used is reduceall, while in the Reid- 
Wittkopf package the procedure is called Clearderivdep. What these procedures do is reduce 
every equation with respect to every other equation until no further reductions are possible. With 
linear systems this results in no loss of information or solutions. The simplification procedure is 
also applied to the output of Kolchin-Ritt (or StandardForm in the Reid-Wittkopf package). Using 
a lexicographic ordering based on ^2 > </> > ?i and x > t > u, (recall in these equations, u is an 
independent variable), the result of simplification, Kolchin-Ritt and further simplification appears 
in the left hand column of Table 1. Clearly, it is necessary to run the calculation again having 
set a + f3 = Q. The result of this second calculation appears in the right hand column of Table 
1. It will be necessary to recalculate the determining equations with the special values of the 
parameters from the beginning if the package used removes factors in the parameters. The output 
is considerably simpler to solve than the original set of equations, while the solutions sets are the 
same. The solutions of these determining equations are given in Table 2. 

We shall now describe symmetry reductions that are generated by these infinitesimals which are 
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a + /3y^0 


a + /3 = 


{(3 + a)it = 

a(/3 + a)(/3</)t + 2e.) = 

Ti - 2e^ = 


/J^xx = 
S^U = 
/J^tt = 

/?</), + 2^4 = 
/3(/?0t + 2ix) = 

Ti - 2C:, = 



Table 1: The classical determining equations ( p.4| ) in Standard Form 





a + /3/0,a/0,/3/0 


a + /3 = 


r 



KiX + K2 

2Kit + K3 

-2Kit/l3 + K4 


KlX + K2t + K3 

2K;it + K4 

-2{Klt + K2x)/P + K5 



Table 2: Classical Infinitesimals 



obtained by solving solving the invariant surface condition ^{x, t, u)ux +t{x, t, u)ut — 0(x, t, u) = 0. 
In the following /ij are constants of integration. 

Case A.l. a + /? / 0. 

(i) If Ki / 0, then we set ki = 1 and ^2 = ^3 = 0, and obtain the symmetry reduction 

u{x,t) = w{z)-t/l3 + \K4\D.t, z = x/t^/^ (3.5) 

where w{z) satisfies 

w"" -\{a + f3)zw'w" - ^a{w'f + (iz^ + ^pK^^ w" + |W - ^^4 = 0. (3.6) 



(ii) If Ki = 0, then we set K3 = 1, and obtain the symmetry reduction 

u{x, t) = w{z) + K^t, Z = X — K2t, 

where w^z) satisfies 

w"" = (a + P)k2w'w" — {1 + K2 + Pk4)w' . 

Setting W = w' , this equation can be integrated twice to yield 

(W'f = |(a + P)k2W^ -{1 + 4+ Pn4)W^ + ^iW + ^lo, 

which is solvable in terms of the Weierstrass elliptic function p{0;g2,g3) (cf. 
Case A.2. a + /? = 0. 



(3.7) 
(3.8) 

(3.9) 
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(i) If Ki / 0, then we set ki = 1 and K3 = K4 = 0, and obtain the symmetry reduction 

U{x,t)=w{z)+ {Kl-l)t-2K2X /P+^K^lnt, 2 = (x + K2t)/t^/^ (3.10) 

where w{z) satisfies 

W"" + ^Piw'f + (iz2 + l^^g) ^" + ^zw' - iK5 = 0. (3.11) 

(ii) If Ki = 0, then set K4 = 1, and we obtain the symmetry reduction 



u{x,t) = w{z) — 2k2 Xt — ^K2t — -^K^t + K^t //3, Z = X — 2K2t — K^t, (3-12) 

where w^z) satisfies the hnear equation 

2k2Z + (2K2K5 — K3 — 1) w" — K2w' — 2K2K3//?. (3.13) 



ini 

w 



3.2.2.B Nonclassical Method. The determining equations for the nonclassical method are 
given in Example 2.4.2. As in the classical case, we need to do the cases a + /3 / and a + /? = 
separately. We assume ^ 7^ 0. 

Case B.l: a + /3 7^ 0. Take a lexicographic ordering based on > ^ and x > t > u. The 
first step of the Reid strategy involves reducing with respect to small linear equations, in this case 
2(l)xu — 3^xx = 0, then calculating the compatibility conditions of each pair of determining equations. 
This leads to several conditions in ^ alone, including 

(2/3 - 3a) (2/3 + a) (3a + 4.l3)iCuC - 2^') = 0. (3.14) 

The precise coefficients in front of this equation depends on the order in which pairs are chosen to 

be cross-differentiated; this affects which equations appear before others and thus which equations 

are used to do the reduction of new compatibility conditions. Continuing with the Reid strategy 

on the equations for ^ yields for generic a and /3 that ^xx = CttS, — 2Q = ^xtC — CxCt = 0. Choosing 

all the special values of a for which ( |3.14] ) is zero and redoing the calculation leads to the same 

conclusion so that 

0-X + K3 
C = —-, , (3.15) 

Kit + K2 

for all values of a and /?, where a is either or 1. 
Subcase B.1.1: ^x / 0. We set a = 1. 



Reducing the determining equations by ( 3.15| ) yields one small equation (p^ = (/i^i— 2) /{Kit + K2) 



and reducing the remaining equations by this still leaves three lengthy equations for (j). In this 
situation, the DirectSearch strategy is beneficial. Using a DirectSearch strategy with the equation 
for (j)^ and each of the lengthy equations leads to two consistency conditions (ki — 2)(/3ki — /3+2a) = 
and a{Ki — 2){2Kia — 4a + (3{k1 — 5ki + 2)) = 0. Needless to say, inserting each subcase in the 
parameters leads to considerable simplification. Continuing in this manner we obtain the solutions 
for each case. There are two cases which yield solutions. 

(i) For all a and /?, with ki = 2 we obtain the infinitesimals 

X + K3 -2t + K3 

these correspond to classical symmetries. 
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(ii) If /? = 2a and ki = 0, then we obtain the infinitesimals 

X + K4 



K2 aK2 

which do not correspond to classical symmetries. 



+ 2 ^ ^ ' +K3, 



K^a 



(3.17) 



Setting K2 = 1/k and K4 = in ( [3. 171) yields the nonclassical reduction 

m(j;, t) = ti;(z) exp(— 2/i:t) + f Kx — tj/(2a) + K3, 2: = xexp(— Kt) 

where w{z) satisfies 



w 



SKazw'w" + Anaww" + ?>na{w') . 



(3.18) 
(3.19) 



Subcase B.1.2: ^^ 



0. 



In this case we set ^ = 1 /(/tit + K2) . Reducing the determining equations with the solution 
for ^ and calculating the compatibility conditions according to the Reid strategy yields ki = 0. If 
a ^ (3 then unless /3 = we obtain only the trivial solution, (j) = fAo) a- constant. 

If a = /?, we can reduce our equations to the simple set 

Setting (p = ^{y), where y = x + nt and = — 12(<& + k^ + 1)//?) the second of the two equations 
can be written as an ode and integrated twice to yield 



d^ 

dy2 



6$^ + Kiy + K2. 



If Ki 7^ then this equation is equivalent the first Painleve equation PI (cf., |^^), whilst if ki = 
then it is solvable in terms of the Weierstrass elliptic functions. Hence we obtain the nonclassical 
reduction 

u{x^t) = v{y) + w{z), y = X + Kt, z = x — Kt, (3.20) 

where v{y) and w{z) satisfy 

+ (1 + K'^)Vyy + 2K(3VyVy 



■"yyyy 



^yy 



-A, 



(3.21) 



and 



w. 



+ (1 + n^)w^^ - 2k(3w^w^^ = A, (3.22) 

respectively, where A is a "separation" constant. Integrating ( 3.21 , 3^) and setting V = Vy W = 
Wz, yields 



Vyy + (1 + k')V + KpV 



-Xy + fii, 



and 



W,, + (1 + K^)W - kPW^ = Xz + /i2, 



(3.23) 



(3.24) 



respectively, where /^i and fi2 are arbitrary constants. If A ^ then these equations are equivalent 
to the first Painleve equation PI, whilst if A = then they are solvable in terms of Weierstrass 
elliptic functions. In particular, if A = ^1 = ^2 = 0, then equations (|3. 23 , 3.24 ) possess the solutions 



V{y) 



3(1 + K^ 

2k/? 



■ sec 



iJl + K^y 



,,^, , 3(1 + ^2) 2 



1 + K? Z 



(3.25) 
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Thus we obtain the exact solution of ( 2.29 ) with a = P given by 



u{x,t) 



Kfi 



tan 



\Jl + K^{x + Kt + 5l] 



H 7^ tan 



\sJl + K'^{x-Kt + 52) ■ 
(3.26) 



Making the transformation x — > ix, t ^ it and u ^ iu in (2.29) with a = (3 yields 

Utt + Uxx + a{UxUxt + UtUxx) - Uxxxx = 0. 

Thus from ( |3.26| ) we obtain the exact solution of ( p. 271) 



(3.27) 



u{x,t) 



3/1+7? 
kP 



tanh 



^Jl + K'^{x + Kt + 5i] 



3yr+7? 

Kf3 



tanh 



Wl + k2 (x - Kt + (52) 



(3.28) 

A plot of this solution is given in Figure 1 and a plot of its derivative with respect to x in Figure 
2. Figure 1 shows that the solution looks like the elastic interaction of a "kink" and an "anti- 
kink" solution. Figure 2 looks like the elastic interaction of two "soliton" solutions. These are of 
particular interest since such solutions are normally associated with integrable equations, whereas 
they arise here for a nonintegrable equation. Furthermore, to our knowledge, this is the first time 
that a "two-soliton" solution has arisen from a nonclassical reduction. Normally such solutions are 
associated with so-called Lie-Backlund transformations (cf., |||]). 

Case B.2: a + /? = 0. The analysis of the determining equations in this case is similar but more 
complicated. We give here only the solution that does not correspond to a classical symmetry. 



^ = Klt + K2, 



u 



+ 



2{koKi — K2)X + "init^ + 2KiK2t^ + 2^2* + K3 



(t + Ko) 

Setting kq we obtain the nonclassical reduction 

u{x,t) = tw{z) + (SK2X + Kit — 6K2^ 



a{t + Ko) 



/q + K3, 



Z = X — i/tlt — K2t, 



where w[z) satisfies 



w"" + a 



{w 



/^2 



WW 



+ (1 + kI)w" - 3kiw' + 2Kl/a = 0. 



(3.29) 



(3.30) 



It is interesting to note that the values of the parameters for which nonclassical reductions were 
found, not corresponding to classical symmetry reductions, were precisely f3 = —a, a and 2a. We 
will not analyse the r = case of the determining equations here. 
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Appendix 1 

We give here the macsyma input files used to calculate the determining equations for nonclassical 
reductions, in the r = 1 case, of the equation ( |2.29| ). They show how to adapt a package designed 
to calculate classical symmetries. 

The batchfile pbql.case containing the macsyma commands to implement the program SYM- 
MGRP.MAX is 

batchloadC'symmgrp.max") ; 

writef ileC'pbq.out") ; 

/* PERTURBED BOUSSINESQ EQUATION - NONCLASSICAL (TAU = 1) */ 

/* u_{tt]- + u_{xx} + a u_x u_{xt]- + b u_t u_{xx]- + u_{xxxx} =0 */ 

batch ("pbql.dat"); 

symmetry (1,0,0) ; 

printeqn(lode) ; 

saveC'lodepbql .Isp" ,lode) ; 

closef ile() ; 



This file in turn batches the file pbql.dat which contains the requisite data about (p.29): 



p:2$ 

q:l$ 

m:l$ 

parameters: [a,b]$ 

warnings : true$ 

sublisteqs: [all]$ 

subst_deriv_of _vi : true$ 

inf o_given : true$ 

highest_derivatives : all$ 

depends ( [etal,eta2,phil] , [x[l] ,x[2] ,u[l]] ) ; 

eta2:l; 

ut:phil-etal*u[l, [1,0]] ; 

uxt:diff (phil,x[l])+diff (phil,u[l])*u[l, [1,0]] - 

diff (etal,x[l])*u[l, [1,0]] 

- diff (etal,u[l])*u[l, [1,0]]**2 - etal*u[l, [2,0]] ; 

utt : diff (phi 1 , x [2] ) +dif f (phil , u [1] ) *ut - diff (etal , x [2] ) *u [1 , [1 , 0] ] 

- diff (etal, u[l])*u[l, [l,0]]*ut - etal*uxt; 
el:utt+u[l, [2,0]]+a*u[l, [1,0]] *uxt+b*ut*u[l, [2,0]]+u[l, [4,0]] ; 
vl:u[l,[4,0]]; 

The important thing to note is that SYMMGRP.MAX recognises that etal and phil represent the 
infinitesimals S, and (j) that are to be determined. We refer the reader to the paper by Champagne, 
Hereman and Winternitz ||l^ for an explanation of the syntax used and the purpose of the various 
commands. 



References 



Algorithms for the Nonclassical Method 24 

[1] R.L. Anderson and N.H. Ibragimov, ^^Lie-Backlund Transformations in Apphcations," SIAM, 
Philadelphia (1979). 

[2] D.J. Arrigo, P. Broadbridge and J.M. Hill, J. Math. Phys., 34 (1993) 4692-4703. 

[3] L. Bachmair and B. Buchberger, ACM SIGSAM Bulletin, 40 (1976) 19-24. 

[4] D. Bayer, "Computational Algebraic Geometry, Part One: Basic Tools", Columbia University, 
Lecture Notes, August 1986. 

[5] D. Bayer and M. Stillman, Duke Math. J., 55 (1987) 321-328. 

[6] D. Bayer and M. Stillman, "Macaulay: A System for computation in algebraic geometry 
and commutative algebra", available for Unix and Macintosh computers (1982-1992). ftp 
128.103.28.10 Macaulay/M3.tar. 

[7] G.W. Bluman and J.D. Cole, J. Math. Mech., 18 (1969) 1025-1042. 

[8] G.W. Bluman and S. Kumei, ^^Symmetries and Differential Equations" Appl. Math. Sci., 81, 
Springer- Verlag, Berlin (1989). 

[9] B. Buchberger, ACM SIGSAM Bulletin, 39 (1976) 19-29. 



[lo; 

[11 

[12 

[13; 

[14 
[15 

[ig; 

[17 

[is; 
[19; 

[20; 

[21 
[22 
[23; 



B. Buchberger, G. Collins and B. Kutzler, Ann. Rev. Comput. Sci., 3 (1988) 85-119. 

G. Carra-Ferro, Proc. AAECG5, Menorca, Spain, (1987) 129-140. 

B. Champagne, W. Hereman and P. Winternitz, Gonip. Phys. Comm., 66 (1991) 319-340. 

B. Char, K. Geddes, G. Gonnet, B. Leong, M. Monagan, and S. Watt, "Mapie V Library 
Reference Manual," Springer Verlag (1991). 

P. A. Clarkson, Chaos, Solitons and Fractals, (1994) to appear. 

P.A. Clarkson and M.D. Kruskal, J. Math. Phys., 30 (1989) 2201-2213. 

P.A. Clarkson and E.L. Mansfield, "Symmetry Reductions and Exact Solutions of a class of 
Nonlinear Heat Equations", Physica D (1993) to appear. 

P.A. Clarkson and E.L. Mansfield, in ^^Applications of Analytic and Geometric Methods to 
Nonlinear Differential Equations," [Ed. P.A. Clarkson], Kluwer, Dordrecht (1993) 375-389. 

P.A. Clarkson and P. Winternitz, Physica, 49D (1991) 257-272. 

D. Cox, J. Little, D. O'Shea, "Ideais, Varieties and Algorithms," Undergraduate Texts in 
Mathematics, Springer Verlag (1992). 

P.G. Estevez, "Nonclassical Symmetries and the Singular Manifold Method for the Burgers' 
Equation", Universidad de Salamanca, (1993). 

W.I. Fushchlich, Ukrain. Mat. Zh., 43 (1991) 1456-1470. 

W.I. Fushchhch and R.Z. Zhdanov, Physics Reports, 172 (1989) 123-174. 

W. Hereman, Euroniath Bull., 2 (1993) to appear. 



Algorithms for the Nonclassical Method 25 



[24 
[25 

[26; 

[27; 
[28; 

[29 

[3o; 

[31 

[32 
[33; 
[34; 
[35 

[36; 

[37; 
[38; 
[39; 
[4o; 

[41 

[42; 

[43; 
[44; 

[45 

[46; 



E.L. Ince, ^^ Ordinary Differential Equations," Dover, New York (1956). 

D. Levi and P. Winternitz, J. Phys. A: Math. Gen., 22 (1989) 2915-2924. 

E. Mansfield, "dif f grob2: A symbolic algebra package for analysing systems of PDE using 
maple", (1993), ftp euclid.exeter.ac.uk, login: anonymous, password: your email address, 
directory: pub/liz. 

E. Mansfield and E. Fackerell, "Differential Grobner Bases", Macquarie Univerisity, preprint 
92/108 (1992). 

E. Mayr and A. Meyer, Adv. Math, 46 (1982) 305-329. 

H. Melenk, H.M. Moller and W. Neun, GROEBNER: A Groebner Basis Package", in ''RE- 
DUCE User's Manual, Version 3.4", A.C. Hearn (1991). 

M.C. Nucci, "Interactive reduce programs for calculating classical, non-classical and Lie- 
Backlund symmetries of differential equations", preprint GT Math:062090-051, School of 
Mathematics, Georgia Institute of Technology, Atlanta (1990). 

P.J. Olver, ''Applications of Lie Groups to Differential Equations," Second Edition, Springer 
Verlag (1993). 

P.J. Olver, "Direct reduction and differential constraints", preprint (1993). 

P.J. Olver and P. Rosenau, Phys. Lett, 114A (1986) 107-112. 

P.J. Olver and P. Rosenau, SIAM J. Appl. Math., 47 (1987) 263-275. 

L.V. Ovsiannikov, "Group Analysis of Differential Equations," [Trans. W.F. Ames], Academic, 
New York (1982). 

E. Pucci, J. Phys. A: Math. Gen., 25 (1992) 2631-2640. 

E. Pucci and G. Saccomandi, J. Math. Anal. Appl, 163 (1992) 588-598. 
G.J. Reid, J. Phys. A: Math. Gen., 23 (1990) L853-L859. 

G.J. Reid, Europ. J. Appl. Math., 2 (1991) 293-318. 

G.J. Reid and A. Wittkopf, A Differential Algebra Package for maple, (1993) ftp 
137 . 82 . 36 . 21 login: anonymous, password: your email address, directory: pub/standardf orm. 

L. Robbiano, Lect. Notes in Comp. ScL, 204 (1985) 513-517. 

F. Schwarz, in "Topics in Soliton Theory and Applications" eds. M.J. Ablowitz, B. Fuchsstein- 
ner and M.D. Kruskal, World Scientific, Singapore, 290-299 (1987). 

W. Trinks, J. Number Theory, 10 (1978) 475-488. 

E.M. Vorob'ev, Acta Appl. Math., 24 (1991) 1-24. 

J. Weiss, M. Tabor and G. Carnevale, J. Math. Phys., 24 (1983) 522-526. 

G.B. Whitham, "Linear and Nonlinear Waves," Wiley, New York (1974). 



Algorithms for the Nonclassical Method 26 

[47] E.E. Whittaker and G.W. Watson, ''Modern Analysis," 4th Edition, C.U.P., Cambridge (1927). 

[48] S. Wolfram, "Matiiematica; A System for Doing Mathematics by Computer," 2nd Edition, 
Addison-Wesley, Redwood-City (1991). 



